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Abstract 

We consider a queue where the server is the Euclidean space, and the customers are ran- 
dom closed sets (RACS) of the Euclidean space. These RACS arrive according to a Poisson 
rain and each of them has a random service time (in the case of hail falling on the Euclidean 
plane, this is the height of the hailstone, whereas the RACS is its footprint). The Euclidean 
space serves customers at speed 1. The service discipline is a hard exclusion rule: no two 
intersecting RACS can be served simultaneously and service is in the First In First Out order: 
only the hailstones in contact with the ground melt at speed 1, whereas the other ones are 
queued; a tagged RACS waits until all RACS arrived before it and intersecting it have fully 
melted before starting its own melting. We give the evolution equations for this queue. We 
prove that it is stable for a sufficiently small arrival intensity, provided the typical diameter of 
the RACS and the typical service time have finite exponential moments. We also discuss the 
percolation properties of the stationary regime of the RACS in the queue. 
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1 Introduction 

Consider a Poisson rain on the d dimensional Euclidean space with intensity A; by Poisson 
rain, we mean a Poisson point process of intensity A in W^^^ which gives the (random) number of 
arrivals in all time-space Borel sets. Each Poisson arrival, say at location x and time t, brings a 
customer with two main characteristics: 
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• A grain C, which is a RACS of ^ centered at the origin. If the RACS is a ball with 
random radius, its center is that of the ball. For more general cases, the center of a RACS 
could be defined as e.g. its gravity center. 

• A random service time a. 

In the most general setting, these two characteristics will be assumed to be marks of the point 
process. In this paper, we will concentrate on the simplest case, which is that of an independent 
marking and independent and identically distributed (i.i.d.) marks: the mark (C, cr) of point (x, t) 
has some given distribution and is independent of everything else. 

The customer arriving at time t and location x with mark (C, a) creates a hailstone, with footprint 
X + C in M'^ and with height a. 

These hailstones do not move: they are to be melted/served by the Euclidean plane at the location 
where they arrive in the FCFS order, respecting some hard exclusion rules: if the footprints of 
two hailstones have a non empty intersection, then the one arriving second has to wait for the end 
of the melting/service of the first to start its melting/service. Once the service of a customer is 
started, it proceeds uninterrupted at speed 1. Once a customer is served/hailstone fully melted, it 
leaves the Euclidean space. 

Notice that the customers being served at any given time form a hard exclusion process as no 
two customers having intersecting footprints are ever served at the same time. For instance, if the 
grains are balls, the footprint balls concurrently served form a hard ball exclusion process. Here 
are a few basic questions on this model: 

• Does there exist any positive A for which this model is (globally) stable? By stability, we 
mean that, for all k and for all bounded Borel set . . . , the vector Ni{t), . . . , Nk{t), 
where Nj{t) denotes the number of RACS which are queued or in service at time t and 
intersect the Borel set Bj, converges in distribution to a finite random vector when t tends 
to infinity. 

• If so, does the stationary regime percolate? By this, we mean that the union of the RACS 
which are queued or in service in a snapshot of the stationary regime has an infinite con- 
nected component. 

The paper is structured as follows. In section |3] we study pure growth models (the ground is 
cold and hailstones do not melt) and show that the heap formed by the customers grows with (at 
most) linear rate with time and that the growth rate tends to zero if the input rate tends to zero. 
We consider models with service (hot ground) in section H] Discrete versions of the problems are 
studied in section [5] 

2 Main Result 

Our main result bears on the construction of the stationary regime of this system. 

As we shall see below (see in particular Equations ([T} and ([T6])). the Poisson Hail model falls in 
the category of infinite dimensional max plus linear systems. This model has nice monotonicity 
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properties (see sections |3]and|4l). However it does not satisfy the separability property of ||2], 
which prevents the use of general sub-additive ergodic theory tools to assess stability, and makes 
the problem interesting. 

Denote by the (random) diameter of the typical RACS (i.e. the maximal distance between its 
points) and by a the service time of that RAC. Assume that the system starts at time t = from 
the empty state and denote by the time to empty the system of all RACS that contain point x 
and that arrive by time t. 

Theorem 1 Assume that the Poisson hail starts at time t = and that the system is empty at that 
time. Assume further that the distributions of the random variables Sf" and a are light-tailed, i.e. 
there is a positive constant c such that Ee'^^'' and Ee'^'^ are finite. Then there exists a positive 
constant Xq ( which depends on d and on the joint distribution of ^ and a ) such that, for any 
A < Ao, the model is globally stable. This means that, for any finite set A in W^, as t ^ oo, the 
distribution of the random field {Wf, x £ A) converges weakly to the stationary one. 



3 Growth Models 

Let <I> be a marked Poisson point process in R'^+^: for all Borel sets B of M*^ and a < b, a. r.v. 

[a, b]) denotes the number of RACS with center located in B that arrive in the time interval 
[a, b]. The marks of this point process are i.i.d. pairs (C„, cr„), where C„ is a RACS of M*^ and an 
is a height (in M+, the positive real line). 

The growth model is best defined by the following equations satisfied by H^, the height at location 
€ M'^ of the heap made of all RACS arrived before time t (i.e. in the (0, t) interval): for all 

t>u>0, 

Ht = K+ I ( < + sup Ry - N-{dv), (1) 
Jlu,t) \ y&C^ ) 

where A'^^ denotes the Poisson point process on M"*" of RACS arrivals intersecting location x\ 

and o"^ (resp. CjJ) the canonical height (resp. RAC) mark process of . That is, if the point 
process A^^ has points Tf , and if one denotes by (erf, Cf ) the mark of point Tf , then o"^ (resp. 
CI) is equal to < (resp. Cf)) on ^Tf , Tf^-^). 

These equations lead to some measurability questions. Below, we will assume that the RACS are 
such that the last supremum actually bears on a subset of Q'^, where Q denotes the set of rational 
numbers, so that these questions do not occur. 

Of course, in order to specify the dynamics, one also needs some initial condition, namely some 
initial field H^, with iT^ G M for all x E W^. 

If one denotes by r^it) the last epoch of A^^' in (— oo, t), then this equation can be rewritten as the 
following recursion: 

Hi =HS+ [ [a: + sup Ry - Hf] N%dv) + + sup Hl^,^. - Hf.^,^ , 
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that is 

Hf = fT^^(j) + sup Hl^,^. lr^(t)>0 + -f^olr^{t)<0- (2) 

These are the forward equations. We will also use the backward equations, which give the heights 
at time for an arrival point process which is the restriction of the Poisson hail to the interval 
[— t, 0] for t > 0. Let Hf denote the height at locations x and time for this point process. 
Assuming that the initial condition is 0, we have 



Mf = cj^.(j) + sup H^_^^,,^. o0^.(i) lr-{t)>-u (3) 
V " """^hi^) ~ J 

with T^(i) the last arrival of the point process in the interval [—t, 0], t > 0, and with the 
time shift on the point processes ifTl . 

Remark 1 Here are a few important remarks on these Poisson hail equations: 

• The last pathwise equations hold for all point processes and all RACS/heights ( although one 
has to specify how to handle ties when RACS with non-empty intersection arrive at the same 
time - we postpone the discussion on this matter to section \5^. 

• These equations can be extended to the case where customers have a more general structure 
than the product of a RACS ofW^ and an interval of the form [0, a]. We will call as profile 
a function s{y, x) : M'^ x M"^ — )• M U {— cxd}, where s{y, x) gives the height at x relative to 
a point y; we will say that point x is constrained by point y in the profile ifs{y, x) ^ —oo. 
The equations for the case where random profiles (rather than product form RACS) arrive 
are 

Hf = sup [Hy^u) +s^^(t)(y,x)) l^x(j)>o + i?olr-(t)<o> (4) 

where (t) is the last date of arrival of before time t, with the point process of 
arrivals of profiles having a point which constrains x. We assume here that this point process 
has a finite intensity. The case of product form RACS considered above is a special case with 
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with the point process of arrivals with RACS intersecting x. 

Here are now some monotonicity properties of these equations: 

1. The representation ^ shows that if we have two marked point processes {N^}x and {N^}x 
such that for all x, C (in the sense that each point of is also a point of N^), and 
if the marks of the common points are unchanged, then Hf < Hf for all t and x whenever 
H^ < H^ for all x. 
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2. Similarly, if we have two marked point processes {N^}x and {N'-^}^ such that for all x, 
jyx ^ j\j-x ("Ijj sense that for all n, the n-th point of A^^ is later than the n-th point of 
N^), and the marks are unchanged, then iJf < Hf for all t and x whenever < Hq for 
all X. 

3. Finally, if the marks of a point process are changed in such a way that C C C and a <a, 
then Hf < Hf for all t and x whenever < Hq for all x. 

These monotonicity properties hold for the backward construction as well. 

They are also easily extended to profiles. For instance, for the last monotonicity property, if 
profiles are changed in such a way that 

s{y,x) <s{y,x), Vx,y, 

then Hf < Hf for all t and x whenever Hq < Hq for all x. 

Below, we use these monotonicity properties to get upper-bounds on the Hf and Mf variables. 

3.1 Discretization of Space 

Consider the lattice Z'^, where Z denotes the set of integers. To each point in a; = {xi, . . . , Xd) G 
we associate the point z{x) = {zi{x), . . . , Zd{x)) € Z'^ with coordinates Zi{x) = [xi\ where 
L J denotes the integer-part. Then, with the RACS A centered at point x G M*^ and having diameter 
^, we associate an auxiliary RACS A centered at point z{x) and being the d-dimensional cube of 
side 2[^J + 2. Since A C A, when replacing the RACS A by the RACS A at each arrival, and 
keeping all other features unchanged, we get from the monotonicity property 3 that for alH G M 
and X G W^, 

Hf<Hf''\ 

with H^ the solution of the discrete state space recursion 

= f f^l^(t) + max H^ , . ) l^.(t)>o + -^olf-(t)<o, ^ G (5) 
with (t) the last epoch of the point process 

in (— cx),t). The last model will be referred to as Model 2. We will denote by R the typical 
half-side of the cubic RACS in this model. These sides are i.i.d. (w.r.t. RACS), and if ^'^ has a 
light-tailed distribution, then has too. 

3.2 Discretization of Time 

The discretization of time is in three steps. 
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Step 1. Model 3 is defined as follows: all RACS centered on z that arrive to Model 2 within time 
interval [n — 1, n), arrive to Model 3 at time instant n — 1. The ties are then solved according to 
the initial continuous time ordering. In view of the monotonicity property 2, Model 3 is an upper 
bound to Model 2. 

Notice that for each n, the arrival process at time n forms a discrete Poisson field of parameter A, 
i.e. the random number of RACS arriving at point z G Z*^ at time n has a Poisson distribution 
with parameter A, and these random variables are i.i.d. in z and n. 

Let {R^ ■, cTni)' — 1)2,..., M,^, be the i.i.d. radii and heights of the cubic RACS arriving at 
point z and time n. Let further M = Mq , = Rq ■, and C7j = cTq j. 

Step 2. Let R^'^"'^ be the maximal half-side of all RACS that arrive at point z and time n in 
Model 2, and i?"*"^ = i?^'"""^. The random variables i?^'"""^ are i.i.d. in z and in n. We adopt the 
convention that i?^'"*"^ = if there is no arrival at this point and this time. If the random variable 
('^ is light-tailed, the distribution of R*^ is also light-tailed, and so is that of {R^'^^'f. Indeed, 

so, for c > 0, 
given Ee^^'' is finite. Let 

z,sum \ ^ z 

i=l 

Then, by Similar arguments, cj*"™ has a light-tailed distribution if a-i do. By monotonicity prop- 
erty 3 (applied to the profile case), when replacing the heap of RACS arriving at {z, n) in Model 
3 by the cube of half-side i?^'"*"^ and of height cTn for all z and n, one again gets an upper 
bound system which will be referred to as Model 4. 

Step 3. The main new feature of the last discrete time Models (3 and 4) is that the RACS that 
arrive at some discrete time on different sites may overlap. Below, we consider the clump made by 
overlapping RACS as a profile and use monotonicity property 3 to get a new upper bound model, 
which will be referred to as Boolean Model 5. 

Consider the following discrete Boolean model, associated with time n. We say that there is a 
"ball" at z at time n if > 1 and that there is no ball at z at this time otherwise. By ball, we 
mean a Lqo ball with center z and radius Rn^""^. By decreasing A, we can make the probabiUty 
p = P(M^ > 1) as small as we wish. 

Let C,f^ be the clump containing point z at time n, which is formally defined as follows: if there is a 
ball at {z, n), or another ball of time n covering z, this clump is the largest union of connected balls 
(these balls are considered as subsets of here) which contains this ball at time n; otherwise, 
the clump is empty. For all sets A of the Euclidean space, let L{A) denote the number of points 
of the lattice Z'^ contained in A. It is known from percolation theory that, for p sufficiently small. 



^sum 0,sum 

and a = (Tq 
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this clump is a.s. finite [6] and, moreover, L{C^) has a light-tailed distribution (since (ij™'*^)'^ is 
light-tailed) [5 |. Recall that the latter means that Eexp(cL(CQ)) < oo, for some c > 0. 

Below, we will denote by Ac the critical value of A below which this clump is a.s. finite and 
light-tailed. 

For each clump C^, let be the total height of all RACS in this clump: 

The convention is again that the last quantity is if C*^ = 0. We conclude also that has a 
light-tailed distribution. 

By using monotonicity property 3 (applied to the profile case), one gets that Boolean Model 5, 
which satisfies the equation 

K = a'n+ max (6) 

with the initial condition Hq = a.s., forms an upper bound to Model 4. Similarly, 

§^ = aii+ max ml_^o9-\ (7) 

S/GCiiU{4 

where is the discrete shift on the sequences {af., C|}. By combining all the bounds constructed 
so far, we get: 

Hf < H^l^'^ and Mf < MjJ^"^^ a.s. (8) 

for all X and t. 

The drawbacks of ^ are twofold: 

(i) for all fixed n, the random variables {C^}z are dependent. This is a major difficulty which 
will be taken care of by building a branching upper-bound in subsections 13.3.11 and 13.3.21 
below. 

(ii) for all given n and z, the random variables and are dependent. We will take care of 
this by building a second upper bound model in subsection 13.3.31 below. 

Each model will bound (O from above and will hence provide an upper bound to the initial con- 
tinuous time, continuous space Poisson hail model. 

3.3 The Branching Upper-bounds 
3.3.1 The Independent Set Version 

Assume that the Boolean Model 5 (considered above) has no infinite clump. Let again be the 
clump containing x G Z"^ at time n. For x 7^ y G Z'^, either = Cn or these two (random) 
sets are disjoint, which shows that these two sets are not independent^ The aim of the following 

'Here "independence of sets" has the probabihstic meaning: two random sets Vi and V2 are independent if P(Vi = 
Ai,\/2 = A2) = P(Vi = yli)P(V2 = ^2),forall4i,42 C Z'^. 
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construction is to show that a certain independent version of these two sets is "larger" (in a sense 
to be made precise below) than their dependent version. 

Below, we call ($7, J", P) the probability space that carries the i.i.d. variables 
from which the random variables {{Cq, ?o)}zgZ'* built. 

Lemma 1 Assume that A < Ac. Let x ^ y be two points in U^. There exists an extension of the 
probability space {Vt, F, P), denoted by (J^, P), which carries another i.i.d. family 

{/ z,sum r,z,max\~i 

and a random pair (Cq , ) built from the latter in the same way as the random variables {{Cq ^a^)} ^^-^d 
are built from {(ctq'''"'", -Ro''""^)}z6Zd. o.nd such that the following properties hold: 

L The inclusion 
holds a.s. 

2. The random pairs (Cg ,5o) (^o>^) independent, i.e. 

P{dS = A^,a^ G B^,d"Q = A2,al G B^) 

= P{CS = ^1, ?o e B,)F{c"q = A2,a^ e B2) 
= P(5o" = ^i,^o €Bi)P(a^ = ^2,£^e52), 

for all sets A\ , Bi and A2 , B^. 

3. The pairs {Q^, £0) {^o^ ^0) ^^^^ same law, i.e. 

= A,al G B)) = P{CyQ = € B), 

for all sets {A, B). 

Proof. We write for short C^ = Cq and = ^0- Consider first the case of balls with a constant 
integer radius R = R"^"-^ (the case with random radii is considered after). Recall that we consider 
Loo -norm balls in W^, i.e. ci-dimensional cubes with side 2R, so a "ball B^ centered at point 
X = {xi, . . . , Xd)" is the closed cube x + [—R, +R]'^. 

We assume that the ball B^ exists at time with probability p = P(M > 1) G (0, 1) indepen- 
dently of all the others. Let = B^ if B^ exists at time and = 0, otherwise, and let 
= I{E^ = B^) be the indicator of the event that B^ exists (we drop the time index to have 
lighter notation). Then the family of r.v.'s {a^j^igz^ i-i-^- 

Recall that the clump C^ , for the input {a^}, is the maximal connected set of balls that contains 
X. This clump is empty if and only if = 0, for all y with doo{x, y) < R. Let L{C^) denote the 
number of lattice points in the clump C^, < L{C^) < 00. Clearly, L{C'^) forms a stationary 
(translation-invariant) sequence. 
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For all sets A c Z'^, let 

Int{A) = {xeA: B"" C A}, and Hit{A) = {xeZ'^: B''nA^$}. 
For A and x,y e A,v/e say that the event 

Int{A),{a^} ^} 

occurs if, for the input {a"}, the random set E"^ = U^G7nt(A) coimected and both x and y 
belong to E^. 

Then the following events are equal: 



zeA V ^' L J z€Hit{A)\Int{A) 

Therefore, the event {C^ = A} belongs to the sigma-algebra generated by the random 

variables {a^,x € Hit{A)}. Let also -^^j^(^) be the sigma-algebra generated by the random 
variables {a^,a^,x G Hit{A)}. 

Recall the notation cjo'^""' = a^j. We will write for short = (To'^"'". Clearly = if 

= 0, and the family of pairs {{a^,a^)} is i.i.d. in 2; G Z''. 
Let {{al, a^)} be another i.i.d. family in 2; G Z'' which does not depend on all random variables 
introduced earlier and whole elements have a common distribution with (q;°,(j°). Let (il, P) 
be the product probability space that carries both {{a^,a^)} and Introduce then a 

third family {(a^,gf )} defined as follows: for any set A containing x, on the event {C^ = A} we 
let 



ia^iA),a^iA)) = 



ial,ai) if zeHit{A) 
(a^,(T^), otherwise. 



When there is no ambiguity, we will use the notation {a^,gf) in place of ia^{A),gf[A)). First, 
we show that {{a^,gf)} is an i.i.d. family. Indeed, for any finite set of distinct points yi, ■ ■ ■ ,yk, 
for any — 1 -valued sequence ii,. . . and for all measurable sets Bi,. . . ,Bk, 

P(a^^ = ij,ay^ G Bj,j = l,...,k) = ^P(C^ = A,ay^ = a* g Bj,j = l,...,k) 

A 

= Yl ^(^"^ = ^' = ^ ^^^yj ^ Hit{A) and a^/i = i-^ay^ G Bj,yj G {Hit{A)y) 

A 

= XI ^l^"" = ^)P("*' = e G i7zt(^))PK^ = ij,ayi G G 

A 

k 

= J2P{C^ = A)llPia'=ij,a'€B,) 
A j=l 

k 

= l[F{a^ = z„a'GBj). 

Notice that the sum over A is a sum over finite A. This keeps the number of terms countable. This 
is Ucit due to assumption on the finiteness of the clumps. 
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Let C_ be the clump of y for {a^} and let = X^^^gf o^^- We now show that the pairs (C^, a^) 

and (C^ , a^) are independent. For all sets A, let J^^ be the sigma-algebra generated by the random 
variables 

and let C^(^) be the clump containing 2/ in the enviroimientQ;"^. Letalsoa^(A) = X^^gc^^^l^)- 
Clearly, (a^"^^ cr^"^)) is also an i.i.d. family. Then, for all sets Ai,Bi and A^, B2, 

P(C^ = G Bi,C' = A2,ay G ^2) 

= P(a" = G Bud' (A,) = A2,ay{Ai) G B2) 

= P(a- = ^1,?^ G B,)P{cnA,) = A2,?n^i) e ^2) 

= P{C'' = Ai,a'' e Bi)P{Cy = A2,ay e B2). (9) 

The second equality follows from the fact that the event {C^ = Ai,a^ G Bi} belongs to the 
sigma-algebra •^^^^(y^^) whereas the event {C_ {Ai) = A2,a^{Ai) G -62} belongs to the sigma- 
algebra T^^, which is independent. The last equaUty follows from the fact that {a^'^^\a^^^^} is 
an i.i.d. family with the same law as {a^,a^}. 

We now prove the first assertion of the lemma. If = , then the inclusion is obvious. 
Otherwise, Q = and if = A, the size and the shape of depend only on {a", u G 
{Hit{A)Y}. Indeed, on these events. 



yint{A%{a^}- 

Then the first assertion follows since, first, the latter relation is determined by {a", u G Int{A'^)} 
and, second, Int{A^) = {Hit{A)Y. We may conclude that C_ {A) ~D because some al,z G 
Hit{A) \ Int(A) may take value 1. 

Finally, the second assertion of the lemma follows from the construction. 
The proof of the deterministic radius case is complete. 

Now we turn to the proof in the case of random radii. Recall that we assume that the radius R of 
a Model 2 RACS is a positive integer- valued r.v. and this is a radius in the L^o norm. For x G Z'^ 
and A; = 1, 2, . . ., let S^'^ be the Loo -norm ball with center r and radius k. Recall that Mq'^ is the 
number of RACS that arrive at time 0, are centered at x and have radius k. Then, in particular, 

Rq = maxjfc : Mq > 1). 

Let a^''' be the indicator of event {Mq ''^ > 1} and E^''' a random set, 

^x,k ^ ^x,k ^x,k ^ ^ ^x,k ^ 0^ otherwise. 

Again, the r.v.'s a^'*^ are mutually independent (now both in x and in k) and also i.i.d. (in x). 
For each A C Z"^, we let Int2{A) = {{x,k) : x e A,k e N,^^'*^ C A} and Hit2{A) = 
{{x, k): x£A,k£n, B'^^^ n ^ 0}- 
For X, y G ^, we say that the event 

Ini2(^),{a«.'} A 
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occurs if, for the input {a^}, the random set = |J^^ k)eint{A) ^^'^ connected and both x 
and y belong to E^. 

Then the following events are equal: 

Therefore, the event {C^ = A} belongs to the sigma-algebra J^Hit2{A) generated by the random 

variables {a^'^, (x, k) S Hit2{A)}. For x £ Z'^ and for = 1, 2, . . ., we let cr^''^ = Ylj=i '^o j 
where the sum of the heights is taken over all RACS that arrive at time 0, are centered at x and have 
radius k. Clearly, the random vectors (a"''^, cr"''^) are independent in all x and k and identically 
distributed in x, for each fixed k. 

Let {(a*''^, (7^''^)} be another independent family of pairs that does not depend on all random 
variables introduced earlier and is such that, for each k and x, the pairs {a^'^ ,a^'^) and {a^'^,a^'^) 
have a common distribution. Let be the product probability space that carries both 

{(a^'^, (7^'^=)} and {(a^''', a^''')}. Introduce then a third family {(a^'^, a^''^)} defined as follows: 
for any set A containing x, on the event {C^ = A} we let 



\aj,a'/) if (zj) eHit2iA) 
(a^'', cr^''), otherwise. 



The rest of the proof is then quite similar to that of the constant radius case: we introduce again 
C_ , which is now the clump of y for {a^''} with the height = J2k Ylzecy show 
that the random pairs (C^, a^) and (C , d_^) are independent and finally establish the first and the 
second assertions of the lemma. 

□ 

We will need the following two remarks on Lemma [T] 

Remark 2 In the proof of Lemma |7] the roles of the points x and y and of the sets and 
are not symmetrical. It is important that is a clump while from V = C^, we only need the 
following monotonicity property: the set V \ is a.s. bigger in the environment {a^} than in the 
environment {a^}. One can note that any finite union of clumps also satisfies this last property. 

Remark 3 From the proof of Lemma\J\ the following properties hold. 

1. On the event where Cq and Cg are disjoint, we have Cq C and ctq'*"'" = ^""^ a.s., 
for all z G Cq, so that CTq < CTq. 

2. On the event where Cq = Cq, we have ag' = CTq. 

Let us deduce from this that, for all constants > a^, for all z € Cg U Cg , there exists a random 
variable r{z) € {x, y} such that z € Cg (with the convention Cg = Cg anda^ = ag j a.s. and 

max ^ (a" + a^) < a''^'^ + Jg^"^ a.s. 

ue{x,y} : z&C^ 
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In case 2 and case 1 with z G Cq, we take r(z) = x and use the fact that > a^. In case I with 
z G Cq, we take r{z) = y and use the fact that Bq < a^. 

As a direct corollary of the last property, the inequality 

max(a'' + ajj , + a^) < max(a'' + CTq , + '^q) 
holds a.s. Here = i^^^Qy CTq 

We are now in a position to fomiulate a more general result: 

Lemma 2 Assume again that A < Ac- Let S be a set of "L^ of cardinality p > 2, say S = 
{xi , . . . ,Xp}. There exists an extension of the initial probability space and random pairs [Cj^ , ct^' ), 
i = 2, . . . ,p defined on this extension which are such that: 

1. The inclusion 

p p 

U ^ U ^o' a.s. (10) 

i=i i=i 

holds with Cq^ = Co^ 



2. For all real valued constants a^^ , a^^ , ■ ■ ■ , such that a^^ = maxi<j<p for all z G 
[J^=i Cq\ there exists a random variable r{z) G {xi, . . . , Xp} such that z G cj ^ a.s. and 

max (a^^- +5o^) < a''(^) a.s. (11) 
je{i,...,p} : z€Co^ 

In particular, the inequality 

max (a^i +5n-' ) < max (a^i (12) 

holds a.s. with = £q\ 

3. The pairs (Cq ^ , CTq ^ ) , {C_q^ , CTq^ ),•••, {C_q^ , CTq'' ) are mutually independent. 

4. The pairs (Cq ' , CTq* ) and {C_q , ag' ), have the same law, for each fixed i = 2, . . . ,p. 

Proof. We proceed by induction on p. Assume the result holds for any set with p points. Then 
consider a set S of cardinality (p + 1) and number its points arbitrarily, S = {xi, . . . , Xp+i). For 
A fixed, consider the event {Cq ^ = A}. On this event, define the same family (a^'', a;^'') as in the 
previous proof and consider the p clumps D_^^ , • • • , 17^^+^ with their heights, say s^^ , . . . , ^v+i 
for this family. By the same reasons as in the proof of Lemma [U (Cq \ Uq ^ ) is independent of 
(^^2 . . . , (^^f+i, s^f+i). By RemarkEl 

P+i p+i 

By the induction step, 

u . . . u D''^+^ Ca,s. Cq' U . . . U Co''+' , 
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with Cg^ , • • • , C_Q^*^ defined as in tlie lemma's statement and then the first, third and fourth asser- 
tions follow. 

We now prove the second assertion, again by induction on p. If p = 2, this is Remark |3] For 
p > 2, we define Li = {p + 1 > j > 1 : Cq' = Cq'} and we consider two cases: 

1. z G In this case let Li = {1, . . . ,p + 1, } \ Li. Since z ^ C^^ for j G Li and since 

= for all j G Li, we get that ([TT]) holds with r{z) = 1 when using the fact that 

a^i = maxi<j<pa^\ 

2. z ^ Cq^. In this case let = {1 < j < p + 1 : j ^ Li,z e C^^}. We can assume 
w.l.g. that this set is non-empty. Then for all j G L^, we have s^^ > o^^j ^ by Lemma[T]and 
Remark |2] So 

max (a^J + CTq^ ) < max {a^^ + s'^j ) a.s. 
jeLi jell 

Now, since the cardinality of is less than or equal to p, we can use the induction as- 
sumption, which shows that when choosing ii G such that a^'i = max^^-j^^ we 
have 

max(a^^' + s^^) < a^'-(-) +^"^"\ 

with r(z) G and with the random variables {a^^} defined as in the lemma's statement, 
but for a_Q^^ which we take equal to s^'i . The proof in concluded in this case too when 
using the fact that the random variable s^'i is mutually independent of the random variables 
({ctq^ }, CTq^ ) and it has the same law as a^'i . 

□ 

3.3.2 Comparison with a Branching Process 

Paths and Heights in Boolean Model 5 Below, we focus on the backward construction associ- 
ated with Boolean Model 5, for which we will need more notation. 

Let denote the set of descendants of level n of x G M"^ in this backward process, defined as 
follows: 

B^ = d^Uix} 
K+i = LJ Cl^^{y}, n>l. 

By construction, B^ is a non-empty set for all x and n. Let denote the cardinality of B^. 

Let denote the set of paths starting from x = xq G Z*^ and of length n in this backward 
process: xq, xi, . . . ,Xn is such a path if xo,xi, . . . , Xn-i is a path of length n — 1 and x„ G 
C^'^^i U {xn-i}- Let vr^ denote the cardinality of H^. Clearly, df^ < vr^ a.s., for all n and x. 

Further, the height of a path /„ = {xq, . . . , x„) is the sum of the heights of all clumps along the 
path: 

n-l 
i=0 



13 



In particular, if the paths In and 4 differ only by the last points Xn € 'S^n+i ^^'^ ^'n ^ ^^-n+v 
then their heights coincide. 

For z € Z*^, let hn^ be the maximal height of all paths of length n that start from x and end at z, 
where the maximum over the empty set is zero. 

Let H^, n > be the maximal height of all paths of length n that start from x. Then H(n) = 
maxx hn ■ 

Paths and Heights in a Branching Process Now we introduce a branching process (also in the 
backward time) that starts from point x = xq at generation 0. Let (V^j, J, z G U^, n > 0, 
i > 1 be a family of mutually independent random pairs such that, for each z, the pair (V^j, sjj J 
has the same distribution as the pair (Cq U {z}, ?q), for all n and i. 
In the branching process defined below, we do not distinguish between points and paths. 

In generation 0, the branching process has one point: IIq" = {(xq}. In generation 1, the points 
of the branching process are = {(xq, xi), xi G VqI}. Here the cardinality of this set is the 
number of points in Vq° and all end coordinates xi differ (but this is not the case for n > 2, in 
general). 

In generation 2, the points of the branching process are 

1X2" = {(xo,Xi,X2), {xQ,Xi) € nf',X2 G V^^{}. 

Here a last coordinate X2 may appear several times, so we introduce a multiplicity function k2'. for 

z G If^, A;| is the number of (xq, xi, X2) G 11^" such that X2 = z. 

Assume the set of all points in generation n is Ilfj" = {(xq, xi, . . . , and A;fj is the multiplicity 
function (for the last coordinate). For each z with k'^ > 0, number arbitrarily all points with last 
coordinate z from I to and let q{xi,X2, ■ ■ ■ ,Xn) denote the number given to point {xq, . . . ,Xn) 
with Xn = z. Then the set of points in generation n + 1 is 

fin+i = {(a^o, • • • ,a;n,a;„+i : (xq, • • • ,x„) G n^°, x„+i G ynl{x^,...,xS^- 

Finally the height of point (xq, . . . , G is defined as 

~ n-l 
/l(xo, . . . , Xn) = ^ sf;^. , 
1=0 

where % = g(xo, . . . ,Xi). 
Coupling of the two Processes 

Lemma 3 Let xq he fixed. Assume that A < Ac. There exists a coupling of Boolean Model 5 and 
of the branching process defined above such that, for all n, for all points z in the set B^o, there 
exists a point (xq, . . . , x„) G such that Xn = z and hn°'^ < h{xo, . . . , x„) a.s. 
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Proof We construct the coupling and prove the properties by induction. For n = 0, 1, the process 
of Boolean Model 5 and the branching process coincide. Assume that the statement of the lemma 
holds up to generation n. For z G ID^o, let = h^'^ . 

Now, conditionally on the values of both processes up to level n inclusive, we perform the fol- 
lowing coupling at level n + we choose with the maximal and we apply Lemma |2] with 
S = B^o, with in place of xi, and with {C^n}z (resp. {C__n}z) in place of {Cq}z (resp. 
{C_q}z)', we then take 

• Kf,i = ^-n U {z} for all z G D^«, z / z*; 

• = ^-n for all z G D^o, z / z*. 

By induction assumption, for all z G D^°, there exists a (xq, . . . , x„) G 11^° such that x„ = z. 
This and Assertion 1 in Lemma|2]show that if n G ^n+v ^^^^ (^Oi • • • i Xn,u) G H^+i, which 
proves the first property. 

By a direct dynamic programming argument, for all u G D^^i 

We get from Assertion 2 in Lemma |2] applied to the set {xi, . . . , Xp} = {z G ID)^°, u G C^„} that 

By the induction assumption, for all z as above, < /i(xo, • • • , x„) a.s. for some (xq, • • • , x„) G 
with x„ = z. Hence for all u as above, there exists a path (xq, . . . , x„, x^+i) G H^^j^ with 
Xn+i = u and such that h'^^^'^ < /i(xo, • • • , x^, x^+i) with (xq, • • • , x„, x^+i) G 11^'^^ and 
x„+i = u. □ 



3.3.3 Independent Heights 

Below, we assume that the light tail assumptions on ^'^ and a are satisfied (see Section |3TI ). 

In the last branching process, the pairs (V^j, J are mutually independent in n, i and z. How- 
ever, for all given n,i and z, the random variables {Vj^^,s^^), are dependent. It follows from 
Proposition [T] in the appendix that one can find random variables {W^ j, J such that 

• For all n, i and z, 1/f ^ C W.^ ^ a.s. 

• The random sets ^ are of the form z + w^^, where the sequence {w^ ^} is i.i.d. in n, i 
and z. 

• The random variable card(Wg has exponential moments. 
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• For all n, i and z, s^.- <t'^^^ a.s. 

• The random variable ^ has exponential moments. 

• The pairs {W^^,t^ J are mutually independent in n, i and z. 

So the branching process built from the {{W^,^, t^i)} variables is an upper bound to the one built 
from the {(^j, J} variables. 

3.4 Upper Bound on the Growth Rate 

The next theorem, which pertains to branching process theory, is not new (see e.g. Pl). We 
nevertheless give a proof for self-containedness. It features a branching process with height (in 
the literature, one also says with age or with height), starting from a single individual, as the one 
defined in Section l3.3.3[ Let v be the typical progeny size, which we assume to be light-tailed. Let 
s be the typical height of a node, which we also assume to be light-tailed. 

Theorem 2 Assume that A < Ac. For n > 0, let h{n) be the maximal height of all descendants of 
generation n in the branching process defined above. There exists a finite and positive constant c 
such that 

lim sup < c a.s. (13) 

n— >cxD n 

Proof. Let (vi, Sj) be the i.i.d. copies of {v, s). Take any positive a. Let D{a) be the event 

Dia) = IjK >a"}, 

n>l 

with dn the number of individuals of generation n in the branching process. For all c > and all 
positive integers k, let Wc^k be the event < c}. Then 

W,,k ^ {W,,knD{a))\jD{a), 

where D{a) is the complement of D{a). From Chemoff 's inequality, we have, for 7 > 

P{D{a)) = P I U {dn+i > d^ < a^ < n} 

\n>0 
n>l \j=l 



n>l 



n>l 



where (/j(7) = Eexp(7?;). First, choose 7 > such that 92(7) < 00. Then, for any integer 
m = 1, 2, . . ., choose am > max(Ew, 2) such that 

qm = V9(7)e-^"'" < ^. 
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So P{D{am)) < 2"" ^ as m ^ oo. 
For any m and any c, 



and 



jlim sup > c| C D{am) M{lim sup — ^ > c} n D{am) 

P ( {lim sup > c} nD{am)) < y'P(n,c,m) , 

\ n^oo JT- / — ' 



where P(n, c, m) = P (^{^ > c} n L>( 
We deduce from the union bound that, for all m, 

/ n \ 

P(n, c, m) < aJ^P Sj > cn . 

The inequality follows from the assumption that u-family and s-family of random variables are 
independent. Hence, by Chemoff 's inequality, 

P(n,c,m)<a;^(^(5)^e-'^'=^ 

where '4^{6) = Ee*^*. Take 5 > such that ^'('5) is finite and then Cm > such that 

hm = a.mi^{6)e-^'''^ < 1. 

Then 



k& 

Hence for all m, 



h(n) 

limsup — < Cmlo^^^y a.s. 

n lb'-' 

Let /i be a random variable taking the value Cm on the event D{am) \ D{am-i)- Then jj, is finite 
a.s. and 

r h{n) ^ 
limsup < /i, a.s. 

n n 

But limsup„ must be a constant (by ergodicity) and then this constant is necessarily finite. 
Indeed, since 

hin) , h(n)oe-^ 
nmsup > nmsup a.s., 

n n n IT' 

and since the shift 9 is ergodic, for each c, the event {lim sup„ < c} has either probability 1 
orO. □ 
Recall that if Ac is the maximal value of intensity A such that Boolean Model 5 has a.s. finite 
clumps, for any A < Ac. 

Corollary 1 Let M{t) = be the height af € Z'^ in the backward Poisson hail growth model 
defined in Under the assumptions of Theorem |2] for all A < Ac, with Ac > the critical 
intensity defined above, there exists a finite constant k(A) such that 

lim sup — ^ = k(A) a.s. (14) 

t^oo t 

with A the intensity of the Poisson rain. 
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Proof. The proof of the fact that the limit is finite is immediate from bound dS]) and Theorem 2. 
The proof that the limit is constant follows from the ergodicity of the underlying model. □ 



Lemma 4 Let a < Ac, where Ac is the critical value defined above. For all A < Ac, 

k(A) = -nia) , (15) 
a 

Proof. A Poisson rain of intensity A on the interval [0, t] can be seen as a Poisson rain of intensity 
a on the time interval [0, At/a]. Hence, with obvious notation 

M(t,A) =El( — ,a) , 



which immediately leads to (fT5l) . □ 

4 Service and Arrivals 

Below, we focus on the equations for the dynamical system with service and arrivals, namely on 
Poisson hail on a hot ground. 

Let denote the residual workload at x and t, namely the time elapsing between t and the 
first epoch when the system is free of all workload arrived before time t and intersecting location 
X G W^. We assume that H§ = 0. Then, with the notation of Section [3l 

= -^ + ^^'(*)+ ^rHt)] lr-W>0- (16) 

V y^c:^(t) J 

We will also consider the Loynes' scheme associated with ([T6]) . namely the random variables 

= Wf o Ot, 

for all X G M'' and t > 0. We have 

= + (t) + sup W^^^.(. o 9r^^t) I Ir^m-t- (17) 

Assume that Wq = Wg = for all x. Using the Loynes-type arguments (see, e.g., Q or ||8]), it 
is easy to show that for all x, Wf is non decreasing in t. Let 

= lim Wf . 

t— )-oo 

By a classical ergodic theory argument, the limit is either finite a.s. or infinite a.s. Therefore, 
for all integers n and all (xi, . . . , Xn) € W^"-, either = oo for all z = 1, . . . , n a.s. or 
< cxo for alH = 1, . . . , n a.s. In the latter case, 

• {^Sol is smallest stationary solution of ([TtT i: 

• (W^'\ . . . , W^" ) converges a.s. to (W;^ , . . . , ) as t tends to oo. 
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Our main result is (with the notation of Corollary [T): 



Theorem 3 If \ < min(Ac, aK{a) then for all x G W^, < oo a.s. 

Proof. For all t > 0, we say that xq is a critical path of length and span t starting from xq in the 
backward growth model {Mf } defined in © if T!:°(t) < -t. The height of this path is H^" = 0. 
For all t > 0, g > 1, we say that xq, xi, . . . , Xg is a critical path of length q and span t starting 
from xo in the backward growth model {Hf } defined in Q if 



with xi G C^XQ, . and T'^°{t) > —t, and if xi, . . . , Xg is a critical path of length q — 1 and span 
t + r^" (t) starting from xi in the backward growth model {H^^^xq o O^xq }. The height of this 
path is EI^°. 

Assume that = oo. Since Wf is a.s. finite for all finite t and all x, there must exist an 



increasing sequence {tk}, with — > oo, such that W^" > > for all k. This in turn 



implies the existence, for all k, of a critical path of length q^ and span tk, say xq, x^, . . . , x^^ of 
height such that 



Then 



for all k and therefore 



tk 



k(A) > liminf — - > 1. 

fc— >oo tk 



Using ([15]) . we get 

But this contradicts the theorem assumptions. □ 



k(A) = —K(a) > 1 a.s. 
a 



Remark 4 Theorem\l\follows from the last theorem and the remarks that precede it. 

Remark 5 We will say that the dynamical system with arrivals and service percolates if there is a 
time for which the directed graph ofRACS present in the system at that time (where directed edges 
between two RACS represent the precedence constraints between them) has an infinite directed 
component. The finiteness of the Loynes variable is equivalent to the non-percolation of this 
dynamical system. 

5 Bernoulli Hail on a Hot Grid 

The aim of this section is to discuss discrete versions of the Poisson Hail model., namely versions 
where the server is the grid rather than the Euclidean space W^. Some specific discrete models 
were already considered in the analysis of the Poisson hail model (see e.g. S ections [3 . 1 1 and [3^ . 
Below, we concentrate on the simplest model, emphasize the main differences with the continuous 
case and give a few examples of explicit bounds and evolution equations. 
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5.1 Models with Bernoulli Arrivals and Constant Services 

The state space is Z. All RACS are pairs of neighbouring points/nodes {i,i + 1}, i G Z with 
service time 1 . In other words, such a RACS requires 1 unit of time for and simultaneous service 
from nodes/servers i and i + 1. For short, a RACS {i, i + 1} will be called "RACS of type i". 
Within each time slot (of size 1), the number of RACS of type i arriving is a Bernoulli-(p) random 
variable. All these variables are mutually independent. If a RACS of type i and a RACS of type 
i + 1 arrive in the same time slot, the FIFO tie is solved at random (with probability 1/2). The 
system is empty at time 0, and RACS start to arrive from time slot (0, 1) on. 

5.1.1 The Growth Model 

(1) The Graph g{l). 

We define a precedence graph ^(1) associated withp = 1 nodes are all {i, n) pairs where i eZis 
a type and n G N = {1, 2, . . .} is a time. There are directed edges between certain nodes, some of 
which are deterministic and some random. These edges represent precedence constraints: an edge 
from {i, n) to (i', n') means that (i, n) ought to be served after (i', n'). Here is the complete list of 
directed edges: 

1. There is either an edge (i + l,n) w.p. 1/2 (exclusive) or an edge (i + l,n) {i,n) 
w.p. 1/2; we call these random edges spatial; 

2. The edges {i, n) — ^ (z — 1, n — 1), (i, n) — ^ {i, n — 1), and {i, n) ^ + 1, n — 1) exist for 
all i and n > 2; we call these random edges time edges. 

Notice that there are at most five directed edges from each node. These edges define directed 
paths: for Xj = {ij,nj), j = 1, . . . , m, the path xi ^> X2 — ^ . . . — > Xm exists if (and only if) all 
edges along this path exist. All paths in this graph are acyclic. If a path exists, its length is the 
number of nodes along the path, i.e. m. 

(2) The Graph a (p). 

We obtain Q{p) from ^(1) by the following thinning: 

1. Each node of ^(1) is colored "white" with probability l—p and "black" with probability p, 
independently of everything else; 

2. If a node is coloured white, then each directed spatial edge from this node is deleted (recall 
that there are at most two such edges); 

3. For n > 2, if a node {i, n) is coloured white, then two time edges (i, n) — > (i — 1, n — 1) 
and (i, n) (i + 1, n — 1) are deleted, and only the "vertical" one, (i, n) {i,n — 1), is 
kept. 

So, the sets of nodes are hence the same in 0(1) and g{p) whereas the set of edges in Q{p) is a 
subset of that in C/(l). Paths in Q{p) are defined as above (a path is made of a sequence of directed 
edges present in g{p)). The graph Q{p) describes the precedence relationship between RACS in 
our basic growth model. 
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The Monotone Property. We have the following monotonicity in p: the smaller p, the thinner 
the graph. In particular, by using the natural coupling, one can make Q{p) C Q {q) for all p < q; 
here inclusion means that the sets of nodes in both graphs are the same and the set of edges of 
Q{p) is included in that of Q{q). 



5.1.2 The Heights and The Maximal Height Function 

We now associate heights to the nodes: the height of a white node is and that of a black one is 1. 
The height of a path is the sum of the heights of the nodes along the path. Clearly, the height of a 
path cannot be bigger than its length. 

For all (i, n), let H}^ = Hn{p) denote the height of the maximal height path among all paths of 
Q{p) which start from node (i,n). By using the natural coupling alluded to above, we get that 
H^Kp) can be made a.s. increasing in p. 

Notice that, for all p < 1, for all n and i, the random variable is finite a.s. To show this, it 
is enough to consider the case p = 1 (thanks to monotonicity) and i = (thanks to translation 
invariance). Let 

*n,n = > 1 ■ ihn)^ {i- 1, n)} 

and, for m = n — 1, n — 2, . . . , 1, let 

tm,n = > *m+i,n + 1 ' {i,m)^{i- l,m)}. 

Similarly, let 

= maxji < —1 : (i, n) — ^ (i + 1, n)} 
and, for m = n — 1, n — 2, . . . , 1, let 

tm,n = max{z < „ - 1 : {i,m)^{i + l,m)}. 

Then all these random variables are finite a.s. (moreover, have finite exponential moments) and 
the following rough estimate holds: 

n 

Hn<Y.{^tn-il,n) + 
i=l 

5.1.3 Time and Space Stationarity 

The driving sequence of RACS is i.i.d. and does not depend on the random ordering of neighbours 
which is again i.i.d., so the model is homogeneous both in time n = 1,2,... and in space i £ Z. 
Then we may extend this relation to non-positive indices of n and then introduce the measure pre- 
serving time-transformation 9 and its iterates 9'^, — oo < m < oo. So H^o9'^ is now representing 
the height of the node {i,n + m) in the model which starts from the empty state at time m. Again, 
due to the space homogeneity, for any fixed n, the distribution of the random variable i?^ does not 
depend on i So, in what follows, we will write for short 

Hn = Hnip) = Hi, 
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when it does not lead to a confusion. 

Definition of function h. We will also consider paths from (0, n) to (0, 1) and we will denote by 
hn = hn{p) the maximal height of all such paths. Clearly, /i„ < a.s. 

5.1.4 Finiteness of the Growtii Rate and Its Continuity at 

Lemma 5 There exists a positive probability po ^ 2/5 such that, for any p < po, 

limsup Hn/n <C{p)< oo, a.s. (18) 

n— >oo 

and 

hri{p) /n ^ ^{p) a.s. and in Li, (19) 
with "fip) and C{p) positive and finite constants, 7(p) < C{p). 

Remark. The sequence {Hn} is neither sub- nor super-additive. 
Lemma 6 For all p, 

lim sup Hn{p) /n < 2j{p) a.s. (20) 

n—^oo 

Lemma 7 Under the foregoing assumptions, 

lim lim sup Hn/ n = a.s. 

PlO n->oo 

Proofs of Lemmas 5-7 are in a similar spirit to that of the main results (Borel-Cantelli lemma, 
branching upper bounds, and also superadditivity), and therefore are omitted. 

5.1.5 Exact Evolution Equations for tlie Growtii Model 

We now describe the exact evolution of the process defined in SJ5TTTT] We adopt here the continuous- 
space interpretation where a RACS of type i is a segment of length 2 centered in z € Z. 

The variable Hn is the height of the last RACS (segment) of type i that arrived among the set with 
time index less than or equal to n (namely with index 1 < /c < n), in the growth model under 
consideration. If {i, n) is black, then Hn is at the same time the height of the maximal height path 
starting from node (i, n) in G{p) and the height of the RACS («, n) in the growth model. If (z, n) 
is white and the last arrival of type i before time n is k, then HJ^ = Hj,. This is depicted in Figure 

m 

If there are no arrivals of type i in this time interval, then Hn = 0. In general, if is the number 
of segments of type i that arrive in [1, n], then Hn > Let Vn be the indicator of the event that 
(i, n) is an arrival (Vn = 1 if it is black and Vn = otherwise). 

Let en*^^ indicate the direction of the edge between (i, n) and (i + 1, n): we write elf'^^ = r if 
the right node has priority, e^n^^ = I if the left node has priority. 
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The following evolution equations hold: if f^+i = 1, then 



h: 



n+l 



V 
V 



Kri + i)i(er' = ^'<+i = i) 
//;;\ + i)i(er''^ = /,^^ = i) 



[Hi V Hi-^ V Ht+^ + 1) 



and if = then = Hl^. Here, for any event A, 1(A) is its indicator function: it equals 

1 if the event occurs and otherwise. 

The evolution equations above may be rewritten as 
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Figure 1 : Top: A realization of the random graph Q (p) . Only the first 6 time-layers are represented. 
A black node at {i, n) represents the arrival of a RACS of type i at time n. Bottom: the associated 
the heap of RACS, with a visualization of the height ff* of each RACS. 



5.1.6 Exact Evolution Equations for the Model with Service 

The system with service can be described as follows: there is an infinite number of servers, each 
of which serves with a unit rate. The servers are located at points 1/2 + i, — oo < i < oo. For 
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each i, RACS {i, n) (or customer (i, n)) is a customer of "type" i that arrives with probability p 
at time n and needs one unit of time for simultaneous service from two servers located at points 
z — 1/2 and i + 1/2. So, at most one customer of each type arrives at each integer time instant. 
If customers of types i and i + 1 arrive at time n, then one makes a decision, that either i arrives 
earUer or i + 1 arrives earUer, at random with equal probabilities, 

P(customer i arrives earUer than customer i + = P(e^*"'"^ = = 1/2- 

Each server serves customers in the order of arrival. A customer leaves the system after the 
completion of its service. As before, we may assume that, for each (z, n), customer {i, n) arrives 
with probability 1, but is "real"("black") with probabiUty p and "virtual"("white") with probabiUty 
\-p. 

Assume that the system is empty at time and that the first customers arrive at time 1. Then, 
for any n = 1, 2, . . ., the quantity := max(T^ — (n — 1), 0) is the residual amount of time 
(starting from time n) which is needed for the last real customer of type i (among customers 
(i, 1), . . . , (i, n)) to receive the service (or equals zero if there is no real customers there). 
Then these random variables satisfy the equations, forn > 1, — cxd < i < cxd, 

TK+i = (W:|l + l)l(e^^+^ = r,t;;+\ = l,<+i = l) 

V (K+i + 1)1(4"''' = h Vn+\ = 1, <+i = 1) 

V ((T^:-l)+ + I«+i = l)) 

V ((K-i-l)+ + l)I(t;;+i = l) 

V i{wt'-iy + i)iK+i = i). 



Since the heights are equal to 1 (and time intervals have length 1), the last two terms in the 
equation may be simplified, for instance, i{W^~^ — 1)+ + = 1) may be replaced by 

wt'iK^, = 1). 

In the case of random heights {cr^}, the random variables {W^} satisfy the recursions 



V + <+i)i(C = I, <-+\ = 1, <+i = 1) 

V ((wri-i)++<+i)i«+i = i) 

V ((K+i-i)++<+i)i«+i = i). 

The following monotonicity property holds: for any n and i, 

^n+l ° ^~""' < ^> a.s. 

Let 

po = sup{p : T{p) < 1}. 

Theorem 4Ifp<po, then, for any i, random variables Wi^ converge weakly to a proper limit. 
Moreover, there exists a stationary random vector {W^, —oo<i<oo} such that, for any finite 
integers iq <0 < ii, the finite-dimensional random vectors 



n 1 n 1 ■ ■ ■ -I n i n ) 
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converge weakly to the vector 

Theorem 5 Ifp < po, then the random variables 

min{i > : = 0} and max{i < : VF* = 0} 

are finite a.s. 

6 Conclusion 

We conclude with a few open questions. The first class of questions pertain to stochastic geometry 
0: 

• How does the RACS exclusion process which is that of the RACS in service at time t in 
steady state compare to other exclusion processes (e.g. Matem, Gibbs)? 

• Assuming that the system is stable, can the undirected graph of RACS present in the steady 
state regime percolate? 

The second class of questions are classical in queueing theory and pertain to existence and prop- 
erties of the stationary regime: 

• In the stable case, does the stationary solution always have a light tail? At the moment, 
we can show this under extra assumptions only. Notice that in spite of the fact that the 
Poisson hail model falls in the category of infinite dimensional max plus linear systems. 
Unfortunately, the techniques developed for analyzing the tails of the stationary regimes of 
finite dimensional max plus linear systems [3 | cannot be applied here. 

• In the stable case, does the Poisson hail equation (fTTl) admit other stationary regimes than 
obtained from {WJ^}^, the minimal stationary regime? 

• For what other service disciplines still respecting the hard exclusion rule like e.g. priorities 
or first/best fit can one also construct a steady state? 

7 Appendix 

Proposition 1 For any pair (X, Y) of random variables with light-tailed marginal distributions, 
there exists a coupling with another pair {(,,f]) of i.i.d. random variables with a common light- 
tailed marginal distribution and such that 

max(X, Y) < min(^, r]) a.s. 
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Proof. Let Fx be the distribution function of X and Fy the distribution function of Y . Let C > 
be such that Ee*^^ and Ee*^^ are finite. Let C = max(0, X, Y). Since e'^'^ < 1 + e"^-^ + e^^, C 
also has a light-tailed distribution, say F. 

LetF(a;) = 1 — F{x), G{x) = F^^^ (x), aadG{x) = 1 — G{x). Let ^ and be i.i.d. with common 
distribution G. Then E'^ is finite, for any c < C/2. 

Finally, a coupUng of X,Y,^, and r} may be built as follows. Let Ui, U2 be two i.i.d. random 

variable having uniform (0,1) distribution. Then let ^ = G~^{Ui), rj = G~^{U2) and C, = 
min(^, 77). Finally, define X and Y given max(X, F) = ^ and conditionally independent of 

□ 
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